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We follow the time sequence of binary elastic collisions in a small collection of hard-core particles. 
Intervals between the collisions are characterized by the numbers of collisions of different pairs in 
a given time. It was shown previously that due to the ergodicity these numbers grow with time as 
a biased random walk. We show that this implies that for a typical trajectory in the phase space 
each particle has "preferences" that are stable during indefinitely long periods of time. During these 
periods the particle collides more with certain particles and less with others. Thus there is a clearly 
distinguishable pattern of collisions of the particle with other particles, as determined by its initial 
position and velocity. The effect holds also for the dilute gas with arbitrary short-range interactions 
allowing for experimental testing. It is the mechanical counterpart to the classical probabilistic 
observation that "in a population of normal coins the majority is necessarily maladjusted" P]. 
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Chaos makes individual trajectories of closed mechan- 
ical systems generally inaccessible to the direct study. 
Correspondingly, the ergodic theory considers instead of 
the evolution of one initial point in the phase space, the 
evolution of finite volumes in the allowed region of the 
phase space (the microcanonical ensemble). The theory 
typically studies the given mechanical system in terms of 
integrals over initial conditions. In particular, these in- 
tegrals define the different-time correlations of functions 
on the phase space (mechanical variables) and the cen- 
tral limit theorem (CLT), which is based on the decay 
of these correlations in time . The latter says roughly 
that for the ergodic systems the mechanical variables be- 
come independent at significantly separated instants of 
time, where the statistical independence is defined by 
the microcanonical ensemble. Then the quantities repre- 
sentable as time integrals of a mechanical variable must 
perform a random walk at large times. It follows that a 
typical individual trajectory of the system must be such 
that on the coarse-grained time-scale the time-averages 
of all its functions are statistically indistinguishable from 
the random walk. Then by a choice of functions one may 
reveal a wealth of physically relevant information on the 
system. The most immediate example is the trajectory 
itself. That can be represented as a time integral of the 
flow in the phase space. Then, assuming the decay of 
the flow's correlations, at large times the trajectory per- 
forms diffusion over the allowed region of the phase space, 
the fact that underlies the whole approach of the ergodic 
theory to the equilibrium statistical mechanics. 

In this Letter we use a different function of the tra- 
jectory to get insight into the details of collisions in the 
dilute gas with arbitrary short-range interactions. The 
main question wc address is: how much the particles col- 
lide with each other within a flnite interval of time t7 To 
make the question meaningful we consider systems with a 
moderate number of particles where each particle collides 
many times with every other particle in the system at a 
reasonable t. We inquire if these collisions are uniform 



so that every particle collides with every other particle 
more or less the same number of times within t. We find 
that typically the numbers of the particle's collisions with 
other particles deviate from uniformity strongly. Typical 
initial conditions for the position and the velocity of the 
particle define the evolution where the particle collides 
with other particles non-equal numbers of times with the 
difference typically growing in time. Thus chaos implies 
a rather counterintuitive property of the system which 
meaning can be made more detailed for hard-core interac- 
tions. In this case the particles bounce off each other in- 
stantaneously. Chaos can be characterized by noting that 
a slight change in the collision against the convex surface 
of the other particle may result in the large change of 
the particle's position in the future. Then our result says 
that sequences of collisions where the particles collide 
with each other more or less the same number of times 
are highly improbable (the initial conditions that corre- 
spond to such trajectories occupy a small fraction of the 
allowed phase space volume). Rather, generically, the 
particle collides more with certain particles and less with 
others during very long intervals of times. These "pref- 
erences" sooner or later are changed to other preferences 
so that eventually all possible preferences will appear, 
however the expectation time to these events is infinite. 

The study assumes that the evolution consists of the 
free motion of the particles interrupted by the binary 
collisions. For hard-core interactions where the collisions 
are instantaneous, this frame holds also in the dense fluid 
phase, to which the results can be extended. We desig- 
nate the number of collisions of the pair of particles i 
and j within time t by Nij{t). These describe times be- 
tween the collisions and on average Nij{t) are equal to 
the ratio of t to the mean free time between the collisions 
of the pair. These numbers were recently represented as 
the time integrals of the mechanical variable [s^] . Using 
facts that hold for the random walk with probability one, 
predictions were made for individual trajectories. 

We use hard balls (disks) as a basic model of study. 
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Many statements of the ergodic theory, that are often 
assumed for the systems studied in statistical mechan- 
ics, are actually proved for this model 0, 01 • The fa- 
mous Boltzmann-Sinai hypothesis states that systems of 
an arbitrary number iV > 2 of elastic hard balls in a 
rf— dimensional box with periodic boundary conditions 
(torus), d > 2, are ergodic in the phase space region 
where the trivial conserved quantities of the system are 
constant [sj . Today this hypothesis can be considered as 
" almost proved" , see e. g. @ and references therein, and 
also 0, SI- The ergodicity does not demand a large N: 
already a system of two disks in a 2— dimensional torus 
is ergodic [9[ . Thus the systems with a moderate number 
of particles considered here can be assumed ergodic. 

The motion of hard balls generates an ordered sequence 
of pairs, say (1,2) (5,6) (7,8)..., meaning that first the 
particles 1 and 2 collided, then 5 and 6, then 7 and 8 
and so on. In particular, this sequence is the basic re- 
sult of the event-driven simulations (EDS) of this sys- 
tem As the list is determined by the Newton law, 
the subsequent pairs in the sequence are not indepen- 
dent. Nevertheless, due to chaos, one expects that well- 
separated pairs in the sequence are independent, so on 
a large time scale the above sequence gets close to the 
sequence of randorn pairs where the next pair is picked 
at random, cf. [ill, [l2| • This expectation was formalized 
in where it was shown that Nij (t) are statistically the 
same as the biased random walk on the coarse-grained 
time-scale. Below we generalize this result to the dilute 
gas with arbitrary short-range interactions. A way to 
define the sequence of non-instantaneous collisions, that 
may overlap in time in this case, is described below. To 
establish the properties of the collisions described above, 
we use the facts of the probability theory [l[ that give 
robust expression to the fact that the behavior of a sin- 
gle realization of the random walk is very different from 
the ensemble average. The main observation is that for 
arbitrarily large number n of the steps of the random 
walk, the most probable fraction of time that the ran- 
dom walk did not change its sign is one. This is in sharp 
contrast to the intuitive expectation that in the limit of 
large n the fraction of time that a long walk is positive 
(negative) is about 1/2. In fact, the fraction 1/2 is the 
least probable, which produces many surprising phenom- 
ena [Ij. The application to the gas seems to be missing 
from the literature. 

The analysis is performed either for dilute gases of par- 
ticles or for the fluid of hard balls, that collide in the 
three-dimensional box with periodic boundary conditions 
(the generalization to hard walls is straightforward) . For 
the gas we neglect triple collisions using the diluteness. 
The basic observation is that Nij{t) can be represented 
as a time-integral of a certain function F{rij, Vij) on the 
phase space, 
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where (t) is the collision rate of the particles i and 



j. Above rij and Vij are the relative coordinate and ve- 
locity of the particles i and j, respectively. We consider 
F{r, v). We assume that fc— th binary collision of parti- 
cles i and j starts at the time , when the two particles 
approach each other to a fixed threshold distance d and 
ends at time tf^ + TcoI, when they cross the threshold 
while separating {tcoI depends on the collision's details). 
For hard balls d is the particles' diameter and TcoI = 0. 
Thus we assume that the interactions have short-range, 
while for the long-range interactions the analysis needs 
refinement, cf. [13]. Using any function x('") that is con- 
stant for r < d, where it equals 1/2, and vanishes in 
a smooth way for r slightly larger than d, we represent 
F{r,v) — —\r ■ v\x'{r)/r, ci. [3]. We have 

- I \r{t) ■ v{t)\x' [r{t)] dt/r{t) = ^ / x' W)] Vmt. 

Performing the integral from times slightly smaller than 
i*^-, where x['"ii(^)] = and < 0, to times slightly 
larger than t^^ + TcoI, where x[?'ij(i)] again equals zero 
with r.ij > 0, we find that the contribution of these times 
to the above integral is unity. Here we assumed that the 
collision occurs so that ry(t) monotonically decreases to 
the minimum and then starts to increase monotonically. 
This assumption is valid for central forces, while appro- 
priate generalizations can be considered in other cases. 
Thus F{r,v) = —\r ■ v\x'{r)/r can be used in Eq. ([T]). 
For our purposes only the existence oi F{r,v) is needed, 
that allows to use the machinery of the ergodic theory to 
study the collision rate (t). In doing so we will assume 
that the non-analyticity in F{r,v) due to the presence 
of |r • t;| factor does not prevent one from using the re- 
sults of the ergodic theory. This assumption is used in 
the kinetic theory to find the collision frequency [l^ and 
it seems sound. We expect a proof can be derived for the 
hard balls, which is left for future work. 

The first conclusion from the ergodic theory is the (ex- 
pected) existence of the average collision rate 

= lim A^„-(t)/t = lim / ^^j{t')dt'/t^ {F{r,v)),{2) 

where we omitted i, j in the last term since v is the same 
for all pairs of particles and it is equal to the inverse aver- 
age time between their collisions. The angular brackets 
designate the equilibrium averages with respect to the 
microcanonical ensemble defined by the trivial conserved 
quantities of the system (for torus the energy and the 
momentum). Thus the statistics is defined by picking 
the initial condition at random in the allowed region of 
the phase space. While the above applies to a single tra- 
jectory, below we address the trajectories statistically. 

It is the deviations of (t) from its mean v determine 
the effect discussed here. The fluctuations produce 

iim ij,mn^ ^ ij.mn 

t^oo t ^ J 

= lim Tij^rnn{t), Tij^rnn{t)= I [(Cij (O)^mn(i)) " dt. 
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One can write ry,„„(<) = {^ij{0)Nmn{t)) - 
(6j(0))(A^mn(i)) so Tij^^ri{t) Hieasures the devia- 
tion of Nmn{t) from the mean vt for evolutions that 
start from the initial conditions where the particles i 
and j collide. We assume that Tij^rnn exist, so ^ij{t) 
effectively has a finite correlation time Tcor < oo, cf. 
0i This assumption is justified in the considered 
three-dimensional case, since one expects the correlations 
of ^ij to decay at large times as where D is the 

space dimension [3] ■ It follows from the CLT that 
at t ^ Tcor the distribution of Nij(t) is Gaussian and 
determined completely by v and Tij^mn- This allows to 
substitute (t) by -I- Qj , where Qj (t) is a white-noise, 

= V + Cy, {Q^3{t)Cmn{t')) = 2r,y- ™„(5(i - t'). (3) 

We note that the use of the CLT above neglects the tails 
in the distribution. Though these tails can be even al- 
gebraic, see e. g. [l^, this is not a limitation for the 
analysis below, that concerns the most probable events. 

The effective Langevin equations ([3]) were introduced 
in Q for the system of hard balls. They produce statisti- 
cally the same process as in Eq. ([T]), on a coarse-grained 
time-scale much larger than Tcor- For hard balls the re- 
sult applies to the fluid-like state of the system, where all 
the particles are well-mixed and all pairs collide (the sit- 
uation where some particles are trapped would produce 
infinite Tcor)- As the original process Nij{t) is discrete 
and it grows in jumps, a more intuitive form of the re- 
sult can be obtained by embedding Eq. (jS]) in a biased 
random walk with parameters defined by v and Tij_rnn- 

The Langevin equations ([3]) produce non-trivial phys- 
ical conclusions [3i;]. The effective correlations and anti- 
correlations between the particles are described by the 
differences Nij^mn{t) = Nij{t) — Nmn{t) that obey the 
equations of the usual Brownian motion 

-^ij,mn — '^ij,mn^ ^ij.mn — Cij Cmn: (4) 
{'^ij,mn(t)'^ij,mn(t )) — 4 [L^j .jj L^j^^^] S(t ^ ), (5) 

The statements that hold for the Brownian motion with 
probability one were applied to derive the deterministic 
predictions about the behavior of the system Q , concern- 
ing the repetitions of the rare events Nij^mnit) = that 
have infinite expectation time. Here we consider the typ- 
ical behavior of Nij^nin{t) providing the bulk of the par- 
ticles' correlations. For a typical evolution of Nij,in{t) 
governed by Eq. Q, possibly after relatively small ini- 
tial fluctuations, Nij^in{t) will perform a long and large 
excursion to, say, positive values (if the excursion is to 
negative values, consider instead Nin^ij[t\ = —Nij^in[t]). 
This excursion means that the particle i will start col- 
liding more with particle j than with particle n. As the 
excursion progresses the difference between the numbers 
of collisions grows. A typical plot, which details will be 
discussed later, is shown in Fig. [T] It shows that the 
difference of the numbers of collisions of the particle 2 
with the particle and particle 3 grows in time steadily 



with small fluctuations. The observer will see that dur- 
ing a relatively long period of time (during which each 
pair collided about 7 x lO"* times), the particle 2 has a 
steady inclination to collide more with particle than 
with particle 3. This looks like there is an effective at- 
traction between the particles 2 and 0. This attraction 
is not due to the force but due to different initial con- 
ditions that generally produce patterns with widely dif- 
ferent Nij{t). This is reinforced by the fact that if one 
waits long enough the random walk -/V20,23(i) will return 
to the origin iV2o,23(i*) — Q aX t — t,, and the effective 
attraction changes to either attraction or repulsion with 
equal probability. This is a non-equilibrium effect per- 
taining to finite times of observation, rather than the 
infinite observation time that corresponds to equilibrium 
in the frame of the ergodic theory. However, the waiting 
time can be arbitrarily large and its probability den- 
sity function (PDF) has a power-law tail tZ^^'^ see e. g. 

The average is infinite, (t,) = oo, and the observer is 
likely to see stable differences in the patterns of particles' 
collisions during the whole time of the experiment. Gen- 
erally one has that the particle i has "preferences" and 
collides more with particle j, than with particle n. The 
preferences hold for long times and can be considered as 
a long-living, quasi-stationary excitation of the system. 

The quantitative description in terms of relative frac- 
tions of the phase space volume of the initial conditions 
that produce given preferences is convenient to formulate 
by embedding Nij^mn{t) in the discrete random walk and 
using the laws closely related to the arcsine law [||. As- 
suming Nij^mnit) changes by unit steps fixes uniquely the 
time-interval At between the jumps from the diffusion co- 
efficient of Wij^rnn{t) iu Eq. ([5]). Then for the sufficiently 
large number 2n of the steps of the random walk, the arc- 
sine law implies that for < a; < I the probability that 
xn time units the walk was positive, and {l—x)n was neg- 
ative, tends to an n— independent limit (2/7r) arcsin-y/x. 
The law implies that Nij,mn{t) keeps its sign in the in- 
terval {t/2,t) with probability 1/2, independently of how 
large t is. Furthermore, the most probable fractions of 
time during which a random walk is positive are either 
close to zero or to one, and not to 1/2. The probability 
that Nij,mn{t) changed its sign k times decreases with k 
and it is most probable that Nij,mn{t) never changed its 
sign for arbitrarily large t. A closely associated result is 
the statement that the probability that the random walk 
changes sign less than Xy/n times tends to a finite limit 
at n — oo. Thus the frequency of the changes decreases 
with n and the duration of intervals of constant sign in- 
creases in length. Moreover, very early and very late sign 
changes are the most probable ones. These features are 

3/9 

closely related to the t» " tail: the return may take ex- 
tremely long times while the random walk performs the 
corresponding large excursion T]. 

The same laws apply to the difference of the total num- 
bers Ni(t) of collisions of different particles, as it is clear 
from the definition Ni{t) — -^iji^)- For a typi- 
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FIG. 1: Shown is the plot of the difference of the number 
of collisions of particle 2 with particle and particle 3. The 
simulation involved 3.5 x 10® events, consisting either of par- 
ticles' collisions between themselves or of passing a wall. Five 
particles were introduced into the box of size 100 x 100 with 
random initial conditions, where in the used units the r. m. 
s. velocity was set to 100. The particles' radius equals 2. 
Periodic boundary conditions were used. Each particle made 
about 7 X 10^ collisions with other particles. 

cal trajectory certain particles collide more, while others 
collide less for long periods of time. To describe different 
preferences of the particles one may use the K{K — l)/2 
dimensional random walk Nij(t) — vt, where K is the 
number of particles in the system. Finally, to show the 
same behavior holds in D = 2 where Tcor = oo we per- 
formed the EDS, see Fig.[TJ Though the difference Nij^mn 
is small relatively to Nij itself, it is systematic, which 
makes it physically relevant. We also verified that similar 
behavior holds for the hard walls' boundary conditions. 

The ergodic theory was used to obtain a rather non- 
trivial information on the behavior of the individual tra- 
jectory of the dilute gas or the fiuid of hard-core particles. 
We showed that the characteristic feature of the chaos - 



the decay of the correlations in time - produces a special 
pattern in the time sequence of the particles' collisions. 
A single realization of the system's evolution has long 
periods of time during which all particles have their own 
preferences, colliding more with certain particles and less 
with others. Thus the initial positions and velocities of 
the particles determine for each one of them a unique 
distinguishable pattern of collisions with other particles. 
This behavior of the mechanical trajectories of the sys- 
tem does not look easily derivable by other means. 

The effect applies to a system with any as long 
as it is ergodic. Its practical observation (numerical or 
experimental) however demands a moderate number of 
particles in view of the large numbers of collisions of the 
same pairs of particles involved [s^. For large the 
system can be divided into cellSj such that the collisions 
occur mostly within the cell |J0 and the results could 
apply locally in space, as it is to be studied. 

The described effect of the preferences is an example of 
the importance of distinguishing the behavior of the sin- 
gle trajectory and ensembled-averaged behavior for non- 
equilibrium systems, cf. [l7j . An important question it 
poses is to obtain the direct explanation from the me- 
chanical point of view. Qualitatively the same physics 
underlies the effect and the chaotic divergence of the tra- 
jectories. Collisions off the convex surfaces of other par- 
ticles are unstable with respect to small changes making 
the trajectories with uniform Nij too special. The sen- 
sitivity to small changes in the initial conditions implies 
that the overall number of collisions' patterns grows fast 
with time, making the preferences - large deviations of 
collisions' patterns from the mean - typical. Similar ef- 
fects are expected in billiards and the Lorentz gas. 

It is feasible to test our predictions for the dilute gas 
with arbitrary short-range interactions experimentally. 
This would give a test of ergodicity for small systems. 

We thank N. Simanyi, N. Chernov, J. R. Dorfman, S. 
Fishman and T. Gotoh for discussions. 
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